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We present a calculation of the next-to-leading order QCD corrections to the partonic cross sections 
contributing to single-inclusive high-pT hadron production in collisions of transversely polarized 
hadrons. We use a recently developed projection technique for treating the phase space integrals in 
the presence of the cos(2$) azimuthal-angular dependence associated with transverse polarization. 
Our phenomenological results show that the double-spin asymmetry A^ T for neutral-pion production 
is expected to be very small for polarized pp scattering at RHIC and could be much larger for the 
proposed experiments with an asymmetric pp collider at the GSI. 
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I. INTRODUCTION 

The partonic structure of spin- 1/2 targets at the leading- 
twist level is characterized entirely by the unpolarized, 
longitudinally polarized, and transversely polarized dis- 
tribution functions /, A/, and Sf, respectively [1]. Of 
these, the "transversity" distributions Sf remain virtu- 
ally unknown. They are defined [1-4] as the differences 
of probabilities for finding a parton of flavor / at scale fi 
and light-cone momentum fraction x with its spin aligned 
(tt) or anti-aligned (11) with that of the transversely po- 
larized nucleoli: 



6f(x, n) = f n (x, /i) - / iT {x, n) 



(1) 



A program of polarized pp collisions is now underway 
at the BNL Relativistic Heavy Ion Collider (RHIC) [5], 
aiming at further unraveling the spin structure of the 
proton. Collisions of transversely polarized protons are 
hoped to give information on transversity through, e.g., 
the measurement of double-spin asymmetries 



Erfo-(TT) - <MU)] 



[dff(TT) + da(U)] 



da 



for various reactions with observed produced high- 
transverse momentum (pt) or invariant mass. The best- 
studied, and perhaps most promising among these, is 
the Drell-Yan process [2, 6, 7], which offers the largest 
spin asymmetries but whose main drawback is the rather 
moderate event rate. Other reactions, such as high-p^ 
prompt-photon, pion, or jet production, are much more 
copious, but suffer from fairly small spin asymmetries 
[3, 8-11], due to large contributions from gluon-gluon 
and quark-gluon scattering only present in the unpolar- 
ized cross section in the denominator of Eq. (2). 



Very recently, it has also been proposed to extract 
transversity from measurements of Att in transversely 
polarized pp collisions at the planned GSI-FAIR facil- 
ity [12-14] near Darmstadt, Germany. For later stages 
of operations, there are plans to have an asymmetric pp 
collider, with moderate proton and antiproton energies of 
3.5 and 15 GeV, respectively. So far, theoretical work has 
focused on the Drell-Yan process. It was found that the 
expected spin asymmetries could be very large, possibly 
reaching several tens of per cents [14-16]. This can be 
readily understood, because for the GSI kinematics only 
partons with rather large momentum fractions scatter off 
each other, and in pp collisions the relevant lowest order 
(LO) process, qq annihilation, will receive large contribu- 
tions from valence quarks, which are expected to carry 
strong polarization. This appears to make the proposed 
measurements at GSI particularly interesting for learning 
about transversity. 

The theoretical framework for GSI kinematics is some- 
what more involved than for RHIC, since perturbative 
all-order resummations of large logarithmic contributions 
to the partonic cross sections are particularly important. 
For the Drell-Yan process, these have been addressed in 
detail recently in [16]. In any case, information from 
RHIC and from GSI experiments will be complementary, 
due to the very different kinematics accessed, and very 
likely both will be needed to gain sufficient knowledge 
about the 6f over a large range in x. 

In this paper, we perform a detailed study of high-py 
single-inclusive pion production in transversely polarized 
pp and pp collisions. In particular, we derive the next- 
to-leading order (NLO) QCD corrections to the relevant 
partonic cross sections. In general, these are indispens- 
able for arriving at a firmer theoretical framework for 
analyzing experimental data in terms of parton densi- 
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ties. For the calculation wc employ a recently developed 
"projection technique" for treating the phase space in- 
tegrals in the presence of the cos(2<j>) azimuthal-angular 
dependence associated with transverse polarization [11]. 

We will apply our analytical results in phenomenolog- 
ical studies for pp and pp scattering at RHIC and GSI 
energies, respectively. Regarding the latter, to our knowl- 
edge our study is the first to propose accessing transver- 
sity via the process pp — > nX. We hope that such mea- 
surements would be possible with the proposed PAX [12] 
and ASSIA [13] experiments, but this reaction should 
also be of great interest for the PANDA Collaboration 
[17]. As we shall see, the spin asymmetry for pp — > irX 
is expected to be much smaller than that for the Drell- 
Yan process, a drawback that may be compensated for 
by the much higher event rates and, therefore, the much 
better statistical accuracy. We will also find, however, 
that the size of the NLO corrections and the scale depen- 
dence are very significant at GSI energies, so that further 
theoretical work will be needed before one can be confi- 
dent that high-pr pion production may be a useful tool 
to learn about transversity. If so, combined information 
from Drell-Yan and from pp — ► ttX (or other produced 



hadrons) could be useful since the processes probe differ- 
ent combinations of the transversity densities. 

In the next section we will very briefly review the nec- 
essary technical framework for the computation of the 
NLO QCD corrections. Details on the projection tech- 
nique, which is a crucial tool for the calculation, can be 
found in Ref. [11] where it was applied to prompt-photon 
production. In Sec. Ill we will present phcnomcnological 
results for RHIC and GSI energies. We conclude in Sec. 
IV. 



II. TECHNICAL FRAMEWORK 

According to the factorization theorem [18], the fully 
differential, transverse-spin dependent, single-inclusive 
cross section for the reaction AB — > irX for the pro- 
duction of a pion (or any other hadron) with transverse 
momentum pr, azimuthal angle $ with respect to the 
initial spin axis, and pseudorapidity ry reads at NLO ac- 
curacy 



y f ^£ /■ 1 - (1 - V)/Ze dv r 1 

Jl-V+VW z c JVW/z c w (l — v ) JV\ 



d 5a pt f L dz, 

dprdrid® nS ab ^ cX Ji-v+vw z t Jvw/z c 



dw 



VW/vz c w 
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d8a 
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"'-' X ^S(1-W) + "'^ dS& ab^cx( S > V > W >») 



dvd$ 



7T 



dvdwdQ 



(3) 



where the sum is over all contributing partonic channels 
ab — ► cX, AB = pp or pp, and with hadron-level variables 



T —II 

V = H W — — c—tv- ' d_\2 

~ S ' S+T 



, S=(P A + P B ) 2 , 



T = (Pa- P,f , U=(P B - P n ) 2 , 



(4) 



in obvious notation of the momenta. The corresponding 
partonic quantities are given by 



t — u 

v = 1 + - , w = — — - 
s s + t 



S=(Pa+ Pbf , 



t = (Po - Pc) 2 , U= (p b - p c ) 2 . (5) 

Neglecting all masses, one has the relations 

s 



x a XbS , t — —T , u — —U 



x b , 



VW 

X a — , Xfo 



i - v 



vwz c ' (1 — v)z c 

The transversity densities in Eq. (3) always refer to 
those for a parent proton even for AB = pp, i.e., we 



use the charge conjugation property 5/f = 5f?. The 
fact that we are observing a specific hadron in the reac- 
tion requires the introduction of additional long-distance 
functions in Eq. (3), the parton-to-pion fragmentation 



functions D* . The dda 



(i) 

ab—>cX 



are the LO (i = 0) and 
NLO (i = 1) contributions in the cross sections for the 
partonic reactions ab — > cX. Finally, fj, collectively de- 
notes the renormalization and factorization scales, which 
we will always take as equal for simplicity. 

We now give a few technical details of the NLO cal- 
culation. Projection on a definite polarization state for 
the initial partons involves the Dirac matrix 75. It is well 
known that in dimensional regularization, which we will 
use to regularize the ultraviolet, infrared, and collinear 
singularities at intermediate stages of the calculation, the 
treatment of 75 is in general a subtle issue. However, 
owing to the chirally odd nature of transversity, in our 
calculation all Dirac traces contain two 75 matrices, and, 
therefore, using the "HVBM scheme" [19] or a naive, to- 
tally anticommuting 75 in n 7^ 4 dimensions must give 
the same results which is also a useful check for the cor- 
rectness of the calculation. 
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The transverse polarization vectors of the initial 
hadrons give rise to a characteristic dependence of the 
cross section on the azimuthal angle $ of the observed 
particle. In the hadronic center-of-mass system (c.m.s.) 
frame, taking the initial hadrons along the ±z axis and 
their spin vectors in ±x direction, the ^-dependence is 
of the form cos(2<I>). Integration over <I> is therefore not 
appropriate. Keeping $ fixed in the NLO calculation 
is however very cumbersome since standard techniques 
developed in the literature for performing NLO phase- 
space integrations rely on the choice of particular refer- 
ence frames different from the one specified above. In [11] 
we developed a general projection method that involves 
integration over all <E>, thereby allowing to keep the ben- 
efits of the standard phase space integration techniques. 
The trick, and the virtue of our method, is to project out 
the dependence of the matrix elements on the spin vec- 
tors in a covariant way, by multiplying with a covariant 
expression for the cos(2<!>) term, and to then carry out 
the complete phase space integrals. 

To be more specific, we note that because of the 
cos(2$) dependence we have the identity 



d 3 Sc 



= cos(2$) [ 
Jo 



d 3 Sc 



7T dpTdrjdQ' 



(7) 

The cos(2$) dependence actually arises through the co- 
variant expression 



F(p c , s a , s b ) 



s 

tu 



tu 

2 (Pc ■ S a ) {p c ■ S b ) + — (S a ■ S b ) 



(8) 

where the Si (i — a, b) are the initial transverse spin vec- 



tors which satisfy s, • p a 



■ Pb = and 



T(p c , s a , s b ) reduces to cos(2$) in the hadronic c.m.s. 
frame. We may, therefore, use F(p c ,s a ,s b )/Tr instead 
of the explicit cos(2<I>)/7r in the "projector" in the inte- 
grand of Eq. (7). For any contributing partonic channel 
we multiply the squared matrix element for transversely 
polarized initial partons, 8\M\ 2 ab ^ cX , by T{p c , s a , s b )/ir. 
The resulting expression may then be integrated over the 
full azimuthal phase space in a covariant way without 
producing a vanishing result, unlike the case of 5|M| 2 
itself; see Ref. [11] for further details. It is crucial here 
that the other observed ("fixed") quantities, the hadron's 
transverse momentum pt and pseudorapidity 77, are de- 
termined entirely by scalar products (p a -p c ) and (p b -p c ), 



independently of the spin vectors s a , b . 

Our method becomes particularly convenient for treat- 
ing the 2—>3 scattering contributions arising at NLO 
where one has an additional phase space integral over the 
second unobserved parton in the final-state. After apply- 
ing the projection method we can perform all phase space 
integrations by employing techniques familiar from the 
corresponding calculations in the unpolarized and lon- 
gitudinally polarized cases [20, 21]. We note that as a 
non-trivial check on our calculation we have also inte- 
grated all squared matrix elements over the spin vectors 
without using any projector at all. This amounts to in- 
tegrating cos(2$) over all < $ < 2ir, and, as expected, 
the final answer is zero. 

The use of dimensional rcgularization is straightfor- 
ward in all this. Ultraviolet poles in the virtual dia- 
grams are removed by the renormalization of the strong 
coupling constant. Infrared singularities cancel in the 
sum between virtual and real-emission diagrams. After 
this cancellation, only collinear poles are left. From the 
factorization theorem it follows that these need to be 
factored into the parton distribution and fragmentation 
functions. This is a standard procedure which we have 
also described in quite some detail in [11, 21]. We use 
the MS scheme throughout. 

After factorization, we arrive at the final result, the fi- 
nite partonic NLO hard scattering cross sections. There 
are all in all five subprocesses that contribute for trans- 
verse polarization: 



qq - 


- qX, 


qq - 


- qX, 


qq - 


- q'X, 


qq - 


- gx, 


qq - 


- gx, 



(9) 



where at NLO in each case X denotes a one- or two- 
parton final state, summed over all possibilities and in- 
tegrated over its phase space. The first four of these 
reactions are present at LO already. The corresponding 
LO transversity cross sections may be found in [8-10, 22]. 
The last subprocess appears for the first time at NLO. 
For each of the five subprocesses, the NLO expression for 
the transversely polarized cross section can be cast into 
the following form: 



dvdwdQ 
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+H ln(l - w) + /ln(l - vw) + J ln(l - v + vw) + K 



In ? 



1 



L 



hi ■±=2- 

1 — VW 



M 



ln(l — v + vw) 



(10) 



where the "plus" -distribution is defined in the usual way 
over the interval [0,1]. All coefficients in Eq. (10) are 
functions of v and w, except those multiplying the distri- 
butions 5(1— w), 1/(1 — w) + , [ln(l — w)/(l — w)] + which 
may be written as functions just of v. Terms with distri- 
butions are present only for the subprocesses that already 
contribute at the Born level. The coefficients are avail- 
able upon request as a FORTRAN code from the authors. 



III. PHENOMENOLOGICAL RESULTS 

We now present some phenomenological results for 
single-inclusive pion production in transversely polarized 
pp collisions at RHIC {VS = 200 and 500 GeV) and 
asymmetric pp collisions at the planned GSI-FAIR facil- 
ity with proton and antiproton energies of 3.5 GeV and 
15 GeV, respectively. 

Since nothing is known experimentally about transver- 
sity so far, we need to model the Sf for our study. Guid- 
ance is provided by the Soffer inequality [23] 



2\Sq{x)\ < q(x) + Aq{x) 



(11) 



which gives bounds for each 5f. As in [10, 11] we utilize 
this inequality by saturating the bound at some low input 
scale no — 0-6 GeV, choosing all signs to be positive, and 
using the NLO (LO) GRV [24] and GRSV ("standard 
scenario") [25] densities q(x,[/,o) and Aq(x,[Ao)> respec- 
tively. For (i > po the transversity densities 8f(x,n) are 
then obtained by evolving them at LO or NLO. We re- 
fer the reader to [7, 10] for more details on our model 
distributions. We note that we will always perform the 
NLO (LO) calculations using NLO (LO) parton distribu- 
tion functions and the two-loop (one-loop) expression for 
a s . We use the pion fragmentation functions of Ref. [26] 
which has both a LO and an NLO set. They provide a 
very good description of the recent RHIC data on unpo- 
larized neutral-pion production [27] . 

Figure 1 shows our estimates for the transversely po- 
larized single-inclusive pion production cross sections at 
LO and NLO for the two different c.m.s. energies at 
RHIC. We have integrated over the range \rj\ < 0.38 
in pseudorapidity, appropriate for measurements with 
the PHENIX detector. Since only half of the pion's az- 
imuthal angle is covered, we integrate over the two quad- 
rants -7r/4 < $ < 7r/4 and 37r/4 < $ < 57r/4, which 

gives ( /.!^/ 4 + / 3 5 J/4 ) cos(2$)d$ = 2. We have also 

varied simultaneously the factorization/renormalization 
scales \x in Eq. (3) within pt < p, < 4p<r; a significant de- 
crease of scale dependence is observed when going from 
LO to NLO. 
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FIG. 1: Transversely polarized single-inclusive neutral-pion pro- 
duction cross sections in LO and NLO for \J~S = 200 and 500 
GeV collisions at RHIC. The LO results have been scaled by a 
factor of 0.01. The shaded bands represent the changes if the 
scale p is varied in the range pr < p. < 4pT- The lower panel 
shows the ratios of the NLO and LO results for both c.m.s. en- 
ergies. 



The lower part of the Figure 1 displays the so-called 
"if-factor" , defined as usual as the ratio of the NLO to 
the LO cross section, for the scale choice fx = 2pr- Ex- 
cept for small pr, where the NLO corrections lead to a 
significant reduction of the cross section, the A"-factor 
turns out to be rather moderate and close to unity. It is 
known that the A-factor for the unpolarized cross section 
is significantly larger than one at RHIC energies, sec, e.g., 
Fig. 4 in Ref. [21] for = 200 GeV, mostly because of 
large corrections found for gluon-initiated partonic chan- 
nels. Therefore, one expects that the double-spin asym- 
metry Att at RHIC will decrease when going from LO to 
NLO. Indeed, as Fig. 2 shows, this is the case. Here we 
used the CTEQ6M (CTEQ6L1) [28] set of unpolarized 
parton distributions to calculate the corresponding NLO 
(LO) unpolarized cross section. We have chosen the scale 
p,= pt which leads to the largest cross sections in Fig. 1. 
We also indicate in Fig. 2 an estimate of the statistical 
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FIG. 2: Upper bounds for the double-transverse spin asymmetry 
Att corresponding to Fig. 1. The "error bars" indicate the 
statistical accuracy that might be achievable at RHIC (see text). 



accuracy that might be achievable at RHIC, based on 

1 



6A TT ~ , (12) 

with beam polarizations Pa,b of 70%, and an integrated 
luminosity C of 320 and 800 pb~ for c.m.s. energies of 
\/S = 200 and 500 GeV, respectively. <7bin denotes the 
unpolarized cross section integrated over the pr-bin for 
which the error is to be determined. Clearly, the statistics 
would be sufficient to measure even asymmetries as small 
as the ones shown in Fig. 2. However, it is likely that the 
systematic error on spin asymmetries at RHIC will not 
be much smaller than 10~ 3 , in which case it would appear 
to be very difficult to access transversity from Att for 
single-inclusive pion production. We stress again that the 
results shown in Fig. 2 are upper bounds, at least within 
the GRV/GRSV framework with its low input scale for 
the evolution. If the bound in Eq. (11) turns out to be 
not saturated at that scale the asymmetries would be 
even smaller. On the other hand, if we used transversity 
densities that saturate Eq. (11) at a higher scale, say 
(tin ~ 1 GeV, the results for Att would be somewhat 
larger. In any case the measurement is very challenging 
at RHIC. 

We now turn to transversely polarized pp collisions 
with = 14.5 GeV at the planned GSI-FAIR facil- 
ity. We first note that for this rather moderate c.m.s. 
energy the pion transverse momentum can at most reach 
7.25 GeV, at mid rapidity. In our study we integrate over 
— 1 < r7i a b < 2.5, where ?7i a b is the pseudorapidity of the 
pion in the laboratory frame. We count positive rapidity 
in the forward direction of the antiproton. For the asym- 
metric collider option we consider here, ?7i a b is related to 
the c.m.s. pseudorapidity 77 via 



??lab = T\ 



2 En 



(13) 



where Ep, E p are the antiproton and proton energies. 
The rapidity interval we use is roughly symmetric in 
c.m.s. pseudorapidities, < 1.75. 

Figure 3 shows our results for the unpolarized (left) 
and transversely polarized (right) cross sections at NLO 
and LO, as functions of p T . For the calculation in the 
unpolarized case we have chosen the GRV [24] parton 
distributions. This choice is motivated by our ansatz for 
the transversity distributions, for which we had also used 
the GRV densities when saturating the Soffer inequality, 
Eq. (11). Unlike at RHIC energies, at y/S = 14.5 GeV 
and pt of several GeV, rather large momentum fractions 
x a ,j, of the partons are probed in Eq. (3), where the polar- 
ized and unpolarized parton densities for a given parton 
type are expected to become similar [29]. It then ap- 
pears most sensible to use the same parton distributions 
in the unpolarized case that we used when modeling the 
transversity densities. In this way we avoid any artifi- 
cial effects in the NLO corrections and Att induced by a 
mismatch in the x — ► 1 behavior of the parton densities 
used in the calculation. 

The shaded bands in the upper panels of Fig. 3 again 
indicate the uncertainties due to scale variation in the 
range px < M < 4pT- One can see that for both, the 
unpolarized and the polarized cross sections, the scale 
dependence does not really improve from LO to NLO. 
This is a characteristic feature in low-order perturbative 
calculations of cross sections for lower fixed-target en- 
ergies, suggesting that corrections beyond NLO are still 
very significant. Indeed, it was recently shown [30] that 
for inclusive-hadron production in the fixed-target regime 
certain double-logarithmic corrections to the partonic 
cross sections are important at each order of perturbation 
theory, and need to be resummed to all orders to achieve 
an adequate theoretical description. Such a resummation 
will be required in particular in the case we are consider- 
ing here and would be very desirable for the future, along 
with a study of power corrections. We emphasize that 
when the proposed measurements of Att will be per- 
formed, it will be crucial to have precise measurements 
also of the unpolarized cross section, in order to test the 
theoretical framework. Only if the theory is sufficiently 
understood will data on Att become useful for determin- 
ing transversity. Similar comparisons of data [27] and 
theoretical calculations for the unpolarized neutral-pion 
cross section at RHIC have shown an excellent agreement 
even down to fairly low pion transverse momenta, which 
has indeed provided much confidence that the calcula- 
tions based on partonic hard-scattering are adequate, so 
that spin asymmetries measured at RHIC determine the 
spin-dependent parton distributions of the proton. 

The lower parts of Fig. 3 display the corresponding K- 
factors at scale \x = 2px- We note that these decrease as 
Pt increases, which is related entirely to the different be- 
havior of the LO and NLO parton distributions at large 
x. Had we chosen the same parton distributions at LO 
and NLO, the if-factors would actually slightly increase 
with pt, as a result of the large double- logarithmic cor- 
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FIG. 3: Unpolarized (left) and transversely polarized (right) single-inclusive neutral-pion production cross sections at LO and NLO 
at the GSI. The LO results have been scaled by a factor of 0.01. The shaded bands represent the theoretical uncertainty if the 
factorization/renormalization scale /i is varied in the range pr < ^ < 4pr. The lower panel shows the ratios of the NLO and LO 
results in each case, using /i = 2pr- 
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FIG. 4: Upper bounds for the transverse double-spin asymmetry 
Att for single-inclusive neutral-pion production in LO and NLO 
at GSI-FAIR. The "error bars" indicate the expected statistical 
accuracy for bins in pr (see text). 



rections mentioned above, which first arise at NLO and 
are known to enhance the cross section. 

Figure 4 shows upper bounds-again within the 
GRV/GRSV framework-for the double-spin asymmetry 
Att for the GSI-FAIR facility. The scale /i is again set 
to px- We also give expectations for the statistical errors 
that may be achievable in experiment. We have calcu- 
lated these using Eq. (12), assuming an integrated lumi- 
nosity of C — 15Qpb~ , and beam polarizations of 30% 
and 50% for the antiprotons and protons, respectively. 



IV. CONCLUSIONS 

We have presented in this paper the complete NLO 
QCD corrections for the partonic hard-scattering cross 
sections relevant for the double-spin asymmetry Att for 
single-inclusive high-p^ pion production in collisions of 
transversely polarized hadrons. This asymmetry could 
be a tool to determine the transversity distributions of 
the nucleon. Our calculation is based on a largely an- 
alytical evaluation of the NLO partonic cross sections, 
and we have used a projection technique for treating the 
characteristic azimuthal-angle dependence introduced by 
the transverse spin vectors. 

In our phenomenological studies we found that the spin 
asymmetry ^4tt is expected to be very small in pp colli- 
sions at RHIC and even decreases when going from LO 
to NLO, due to a larger A"-factor in the unpolarized case. 
We have also studied Att for possibly forthcoming trans- 
versely polarized pp collisions in an asymmetric collider 
mode at the GSI-FAIR facility. Here, the spin asym- 
metry may be much larger, but it will be crucial in the 
future to investigate the effects of all-order resummations 
of large Sudakov logarithms. Detailed measurements of 
the unpolarized cross sections will be essential for test- 
ing the applicability of the theoretical framework at the 
moderate c.m.s. energies available at GSI-FAIR. 
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